The Late Embryogenesis Abundant Proteins in Soybean: Identification, Expression Analysis, and the Roles of GmLEA4_19 in Drought Stress.

Late embryogenesis abundant (LEA) proteins play important roles in regulating plant growth and responses to various abiotic stresses. In this research, a genome-wide survey was conducted to recognize the LEA genes in Glycine max. A total of 74 GmLEA was identified and classified into nine subfamilies based on their conserved domains and the phylogenetic analysis. Subcellular localization, the duplication of genes, gene structure, the conserved motif, and the prediction of cis-regulatory elements and tissue expression pattern were then conducted to characterize GmLEAs. The expression profile analysis indicated that the expression of several GmLEAs was a response to drought and salt stress. The co-expression-based gene network analysis suggested that soybean LEA proteins may exert regulatory effects through the metabolic pathways. We further explored GnLEA4_19 function in Arabidopsis and the results suggests that overexpressed GmLEA4_19 in Arabidopsis increased plant height under mild or serious drought stress. Moreover, the overexpressed GmLEA4_19 soybean also showed a drought tolerance phenotype. These results indicated that GmLEA4_19 plays an important role in the tolerance to drought and will contribute to the development of the soybean transgenic with enhanced drought tolerance and better yield. Taken together, this study provided insight for better understanding the biological roles of LEA genes in soybean.

Characterization of the soybean KRP gene family reveals a key role for GmKRP2a in root development.

Kip-related proteins (KRPs), as inhibitory proteins of cyclin-dependent kinases, are involved in the growth and development of plants by regulating the activity of the CYC-CDK complex to control cell cycle progression. The KRP gene family has been identified in several plants, and several KRP proteins from Arabidopsis thaliana have been functionally characterized. However, there is little research on KRP genes in soybean, which is an economically important crop. In this study, we identified nine GmKRP genes in the Glycine max genome using HMM modeling and BLASTP searches. Protein subcellular localization and conserved motif analysis showed soybean KRP proteins located in the nucleus, and the C-terminal protein sequence was highly conserved. By investigating the expression patterns in various tissues, we found that all GmKRPs exhibited transcript abundance, while several showed tissue-specific expression patterns. By analyzing the promoter region, we found that light, low temperature, an anaerobic environment, and hormones-related cis-elements were abundant. In addition, we performed a co-expression analysis of the GmKRP gene family, followed by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) set enrichment analysis. The co-expressing genes were mainly involved in RNA synthesis and modification and energy metabolism. Furthermore, the GmKRP2a gene, a member of the soybean KRP family, was cloned for further functional analysis. GmKRP2a is located in the nucleus and participates in root development by regulating cell cycle progression. RNA-seq results indicated that GmKRP2a is involved in cell cycle regulation through ribosome regulation, cell expansion, hormone response, stress response, and plant pathogen response pathways. To our knowledge, this is the first study to identify and characterize the KRP gene family in soybean.

Wavefront aberrations caused by biomechanical effects after Small Incision Lenticule Extraction (SMILE) based on finite element analysis.

To examine wavefront aberrations induced by biomechanical effects after Small Incision Lenticule Extraction (SMILE) surgery. The three-dimensional (3D) finite element models of the human eye were established. By loading the intraocular pressure (IOP), the displacement of the anterior and posterior surface of the cornea was calculated. Then the displacement was converted into the wavefront aberrations by wave-surface fitting. The results showed that the induced wavefront aberrations were noticeable from biomechanical effects after SMILE surgery. The induced higher-order aberrations from the anterior corneal surface included spherical aberration, y-Trefoil, and x-Tetrafoil. Spherical aberration was positively correlated with corrected diopter (D), but x-Tetrafoil and y-Trefoil remained stable. The induced wavefront aberrations from the posterior corneal surface were smaller than those from the anterior corneal surface, and some of the aberrations compensated each other. With IOP increased, defocus and x-Tetrafoil from the anterior corneal surface increased, while y-Trefoil and spherical aberration decreased. The IOP only affected defocus from the posterior corneal surface. In addition, the incision size also had a distinct impact on primary x-astigmatism and x-Trefoil from the anterior corneal surface, and it had a smaller effect on the aberrations from the posterior corneal surface. Therefore, the biomechanical effects increased residual wavefront aberrations after SMILE refractive surgery.

Biomechanical responses of the cornea after small incision lenticule extraction (SMILE) refractive surgery based on a finite element model of the human eye.

PURPOSE: To investigate the biomechanical responses of the human cornea after small incision lenticule extraction (SMILE) procedures, especially their effects of SMILE surgery on stress and strain. METHODS: Based on finite element analysis, a three-dimensional (3D) model of the human eye was established to simulate SMILE refractive surgery procedures. Stress and strain values were calculated by inputting the intraocular pressure (IOP). RESULTS: After SMILE refractive surgery procedures, the stress and strain of the anterior and posterior corneal surfaces were significantly increased. The equivalent stress and strain on the anterior and posterior corneal surfaces increased with increasing diopter and were concentrated in the central area, whereas the values of stress and strain at the incision site on the anterior surface of the cornea were approximately 0. Compared with the anterior corneal surface, the stress and strain of the posterior surface were larger. Increasing IOP caused an approximately linear change in stress and a nonlinear increase in corneal strain. In addition, we found that the incision sizes and direction had less of an influence on stress and strain. In summary, SMILE surgery increased the equivalent stress and strain on the human cornea. CONCLUSIONS: The equivalent stress and strain of the anterior and posterior human corneal surfaces increased after SMILE refractive surgery; these increases were particularly noticeable on the posterior surface of the cornea.

Genome-wide analysis of the soybean root transcriptome reveals the impact of nitrate on alternative splicing.

In plants, nitrate acts not only as a signaling molecule that affects plant development but also as a nutrient. The development of plant roots, which directly absorb nutrients, is greatly affected by nitrate supply. Alternative gene splicing plays a crucial role in the plant stress response by increasing transcriptome diversity. The effects of nitrate supply on alternative splicing (AS), however, have not been investigated in soybean roots. We used high-quality high-throughput RNA-sequencing data to investigate genome-wide AS events in soybean roots in response to various levels of nitrate supply. In total, we identified 355 nitrate-responsive AS events between optimal and high nitrate levels (NH), 335 nitrate-responsive AS events between optimal and low nitrate levels (NL), and 588 nitrate-responsive AS events between low and high nitrate levels (NLH). RI and A3SS were the most common AS types; in particular, they accounted for 67% of all AS events under all conditions. This increased complex and diversity of AS events regulation might be associated with the soybean response to nitrate. Functional ontology enrichment analysis suggested that the differentially splicing genes were associated with several pathways, including spliceosome, base excision repair, mRNA surveillance pathway and so on. Finally, we validated several AS events using reverse transcription-polymerase chain reaction to confirm our RNA-seq results. In summary, we characterized the features and patterns of genome-wide AS in the soybean root exposed to different nitrate levels, and our results revealed that AS is an important mechanism of nitrate-response regulation in the soybean root.

Simulated biomechanical effect of aspheric transition zone ablation profiles after conventional hyperopia refractive surgery.

We studied the effects of the aspheric transition zone on the optical wavefront aberrations, corneal surface displacement, and stress induced by the biomechanical properties of the cornea after conventional laser in situ keratomileusis (LASIK) refractive surgery. The findings in this study can help improve visual quality after refractive surgery. Hyperopia correction in 1-5D was simulated using five types of aspheric transition zones with finite element modeling. The algorithm for the simulations was designed according to the optical path difference. Wavefront aberrations were calculated from the displacements on the anterior and posterior corneal surfaces. The vertex displacements and stress on the corneal surface were also evaluated. The results showed that the aspheric transition zone has an effect on the postoperative visual quality. The main wavefront aberrations on the anterior corneal surface are defocus, y-primary astigmatism, x-coma, and spherical aberrations. The wavefront aberrations on the corneal posterior surface were relatively small and vertex displacements on the posterior corneal surface were not significantly affected by the aspheric transition zone. Stress analysis revealed that the stress on the cutting edge of the anterior corneal surface decreased with the number of aspheric transition zone increased, and profile #1 resulted in the maximum stress. The stress on the posterior surface of the cornea was more concentrated in the central region and was less than that on the anterior corneal surface overall. The results showed that the aspheric transition zone has an effect on postoperative aberrations, but wavefront aberrations cannot be eliminated. In addition, the aspheric transition zone influences the postoperative biomechanical properties of the cornea, which significantly affect the postoperative visual quality.

Genome-wide transcriptional profiling for elucidating the effects of brassinosteroids on Glycine max during early vegetative development.

Soybean is a widely grown grain legume and one of the most important economic crop species. Brassinosteroids play a crucial role in plant vegetative growth and reproductive development. However, it remains unclear how BRs regulate the developmental processes in soybean, and the molecular mechanism underlying soybean early development is largely unexplored. In this study, we first characterized how soybean early vegetative growth was specifically regulated by the BR biosynthesis inhibitor propiconazole; this characterization included shortened root and shoot lengths, reduced leaf area, and decreased chlorophyll content. In addition, the growth inhibition induced by Pcz could be rescued by exogenous brassinolide application. The RNA-seq technique was employed to investigate the BR regulatory networks during soybean early vegetative development. Identification and analysis of differentially expressed genes indicated that BRs orchestrate a wide range of cellular activities and biological processes in soybean under various BR concentrations. The regulatory networks between BRs and multiple hormones or stress-related pathways were investigated. The results provide a comprehensive view of the physiological functions of BRs and new insights into the molecular mechanisms at the transcriptional level of BR regulation of soybean early development.

Identification and characterization of a stachyose synthase gene controlling reduced stachyose content in soybean.

KEY MESSAGE: We identified and characterized a mutant of soybean stachyose synthase gene controlling reduced stachyose content which benefit the soybean seed composition breeding program in the future. It has been shown that in soybean, increased sucrose and reduced raffinose family oligosaccharides would have a positive impact on the world's feed industry by improving digestibility and feed efficiency. We searched for new sources of modified oligosaccharide content in a subset of the USDA Soybean Germplasm Collection and then identified plant introduction (PI) 603176A as having ultra-low stachyose content (0.5%). We identified a 33-bp deletion mutant in the putative stachyose synthase gene (STS gene, Glyma19g40550) of PI 603176A. A co-dominate indel marker was successfully developed from this 33-bp deletion area and was genetically mapped into two F 2:3 populations and a F 4:5 population, which associated with low stachyose content in the progeny lines. These observations provided strong evidence that the STS gene is responsible for stachyose biosynthesis in the soybean plant. Expression of the sts gene remained at the normal level, suggesting the loss of function in the gene is due to defective protein function. This gene-based perfect genetic marker for low stachyose content can be useful for marker-assisted selection in soybean molecular breeding programs.

Molecular cloning, expression profiling and functional analyses of a cDNA encoding isopentenyl diphosphate isomerase from Gossypium barbadense.

Gossypol, a type of plant defence sesquiterpenoid phytoalexin, is synthesized from the MEP (2C-methyl-D-erythritol 4-phosphate) and MVA (mevalonate) pathway in the isoprenoid biosynthetic system. The key step is the isomerization of IPP (isopentenyl diphosphate) to DMAPP (dimethylallyl diphosphate), which is catalysed by IPI (IPP isomerase; EC 5.3.3.2). A full-length cDNA encoding IPI (designated GbIPI) was cloned from Gossypium barbadense by RACE (rapid amplification of cDNA ends). The full-length cDNA of GbIPI was 1205 bp and contained a 906 bp ORF (open reading frame) encoding a protein of 302 amino acids, with a predicted molecular mass of 34.39 kDa and an isoelectric point of 6.07. Amino acid sequence analysis revealed that the GbIPI has a high level of similarity to other IPIs. Southern-blot analysis revealed that GbIPI belongs to a small gene family. Expression analysis indicated that GbIPI expression is highest in stems, followed by leaves, and is lowest in roots, and that the expression of GbIPI could be induced by Verticillium dahliae Kleb, MeJA (methyl jasmonate) and SA (salicylic acid). The functional colour assay indicated that GbIPI could accelerate the accumulation of beta-carotene in Escherichia coli transformants. The cloning and functional analysis of GbIPI will be useful in increasing understanding of the role of IPI in isoprenoid biosynthesis at the molecular level.

Molecular cloning and functional analysis of the gene encoding 3-hydroxy-3-methylglutaryl coenzyme A reductase from hazel (Corylus avellana L. Gasaway).

The enzyme 3-hydroxy-3-methylglutaryl-CoA reductase (HMGR; EC1.1.1.34) catalyzes the first committed step of isoprenoids biosynthesis in MVA pathway. Here we report for the first time the cloning and characterization of a full-length cDNA encoding HMGR (designated as CgHMGR, GenBank accession number EF206343) from hazel (Corylus avellana L. Gasaway), a taxol-producing plant species. The full-length cDNA of CgHMGR was 2064 bp containing a 1704-bp ORF encoding 567 amino acids. Bioinformatic analyses revealed that the deduced CgHMGR had extensive homology with other plant HMGRs and contained two transmembrane domains and a catalytic domain. The predicted 3-D model of CgHMGR had a typical spatial structure of HMGRs. Southern blot analysis indicated that CgHMGR belonged to a small gene family. Expression analysis revealed that CgHMGR expressed high in roots, and low in leaves and stems, and the expression of CgHMGR could be up-regulated by methyl jasmonate (MeJA). The functional color assay in Escherichia coli showed that CgHMGR could accelerate the biosynthesis of beta-carotene, indicating that CgHMGR encoded a functional protein. The cloning, characterization and functional analysis of CgHMGR gene will enable us to further understand the role of CgHMGR involved in taxol biosynthetic pathway in C. avellana at molecular level.

Transformation of taxol-producing endophytic fungi by restriction enzyme-mediated integration (REMI).

The REMI method was used to introduce the plasmid pV2 harboring the hygromycin B phosphotransferase (hph) gene controlled by the Aspergillus nidulans trpC promoter and the trpC terminator into a taxol-producing endophytic fungus BT2. REMI transformation yielded stable transformants capable of continuing to grow on PDA medium containing 125 mug mL(-1) hygromycin B. The transformation efficiency was about 5-6 transformants mug(-1) plasmid DNA. The presence of hph gene in transformants was confirmed by PCR and Southern blot analyses. To the authors' knowledge, this is the first report on the transformation of taxol-producing endophytic fungi by the REMI technique. This study provides an effective approach for improving taxol production of endophytic fungi by the genetic engineering of taxol biosynthetic pathway genes in the future.

Molecular cloning and heterologous expression of a 10-deacetylbaccatin III-10-O-acetyl transferase cDNA from Taxus x media.

A full-length cDNA encoding 10-deacetylbaccatin III-10-O-acetyl transferase (designated as TmDBAT), which catalyzes the acetylation of the C-10 hydroxyl group of the advanced metabolite 10-deacetylbaccatin III (10-DAB) to yield baccatin III, the immediate diterpenoid precursor of Taxol, was isolated from Taxus x media. Heterologous expression of TmDBAT in E. coli demonstrated that TmDBAT was a functional gene. Tissue expression pattern analysis revealed that TmDBAT expressed strongly in leaves, weak in stems and no expression could be detected in fruits, implying that TmDBAT was tissue-specific. Expression profiling analysis of TmDBAT under different elicitor treatments including silver nitrate, ammonium ceric sulphate and methyl jasmonate indicated that TmDBAT was an elicitor-responsive gene. Southern blot analysis suggested that TmDBAT belonged to a small multigene family.

Molecular cloning and expression profile analysis of Ginkgo biloba DXS gene encoding 1-deoxy-D-xylulose 5-phosphate synthase, the first committed enzyme of the 2-C-methyl-D-erythritol 4-phosphate pathway.

Plant diterpenes such as ginkgolides are biosynthesized via the recently discovered 2-C-methyl-D-erythritol 4-phosphate (MEP) pathway. The initial step of the MEP pathway is the formation of 1-deoxy-D-xylulose 5-phosphate (DXP) catalyzed by 1-deoxy-D-xylulose 5-phosphate synthase (DXS, EC: 4.1.3.37), which may thus be considered the first committed step of the MEP pathway for ginkgolides biosynthesis. The full-length cDNA of DXS was isolated and characterized from the gymnosperm plant species, Ginkgo biloba. The full-length cDNA of GbDXS was 2795 bp containing a 2154 bp open reading frame (ORF) encoding 717 amino acids. Comparative and bioinformatic analyses revealed that GbDXS has extensive homology with DXSs from other plant species and, like these, contains a conserved transit peptide for plastid import, histidine residue, a putative thiamine diphosphate-binding site and a transketolase motif. Phylogenetic analysis indicates that GbDXS belongs to the plant DXS1 cluster and suggests it to be more ancient than other plant DXSs. GbDXS was found to be expressed in all tested tissues including roots, stems, leaves, pericarps and seeds. Expression profiling analyses revealed that GbDXS expression was induced by exogenous elicitors including methyl jasmonate, arachidonic acid, acetylsalicylic acid and ceric ammonium sulfate, and showed that the transcription levels were correlated with ginkgolide accumulation, suggesting that DXS might play a regulatory role in ginkgolide biosynthesis in cell culture of G. biloba at the transcriptional level.

[Isolation and characterization of a 2C-methyl-D-erythritol 2,4-cyclodiphosphate synthase gene from Taxus media].

2C-methyl-D-erythritol 2,4-cyclodiphosphate (MEC) synthase (MECS, EC: 4.6.1.12) is the fifth enzyme of the nonmevalonate terpenoid pathway for isopentenyl diphosphate biosynthesis and further Taxol biosynthesis. The full-length MECS cDNA sequence (GenBank accession number DQ286391) was cloned and characterized for the first time from Taxus media, using Rapid Amplification of cDNA Ends (RACE) technique. The full-length cDNA of Tmmecs was 1081 bp containing a 741 bp open reading frame (ORF) encoding a peptide of 247 amino acids with a calculated molecular mass of 26.1 kDa and an isoelectric point of 8.97. Comparative and bioinformatic analyses revealed that TmMECS had extensive homology with MECSs from other plant species. Phylogenetic analysis indicated that TmMECS was more ancient than other plant MECSs. Southern blot analysis revealed that Tmmecs belonged to a small gene family. Tissue expression pattern analysis indicated that Tmmecs expressed constitutively in all tissues including roots, stems and leaves. The cloning and characterization of Tmmecs will be helpful to understand more about the role of MECS involved in the Taxol biosynthesis at the molecular level.

The development of biotechnology education in China*.

From the middle of the 20th century, Chinese scientists have been actively involved in biotechnology. However, biotechnology education in China is a relatively recent phenomenon. This subject has not been addressed at the undergraduate level in a serious way until recently. In the last decade, biotechnology education developed rapidly and reached a new level in Chinese universities. The Chinese scientific establishment is very much aware of the importance of biotechnology and has identified this subject as one of the priority areas. Some universities are taking positive steps toward enhancing biotechnology education. This article focuses on the emergence, as well as the problems and prospects, of biotechnology education in China.

Characterization and expression profile analysis of a new cDNA encoding taxadiene synthase from Taxus media.

A full-length cDNA encoding taxadiene synthase (designated as TmTXS), which catalyzes the first committed step in the Taxol biosynthetic pathway, was isolated from young leaves of Taxus media by rapid amplification of cDNA ends (RACE). The full-length cDNA of TmTXS had a 2586 bp open reading frame (ORF) encoding a protein of 862 amino acid residues. The deduced protein had isoelectric point (pI) of 5.32 and a calculated molecular weight of about 98 kDa, similar to previously cloned diterpene cyclases from other Taxus species such as T. brevifolia and T. chinenisis. Sequence comparison analysis showed that TmTXS had high similarity with other members of terpene synthase family of plant origin. Tissue expression pattern analysis revealed that TmTXS expressed strongly in leaves, weak in stems and no expression could be detected in fruits. This is the first report on the mRNA expression profile of genes encoding key enzymes involved in Taxol biosynthetic pathway in different tissues of Taxus plants. Phylogenetic tree analysis showed that TmTXS had closest relationship with taxadiene synthase from T. baccata followed by those from T. chinenisis and T. brevifolia. Expression profiles revealed by RT-PCR under different chemical elicitor treatments such as methyl jasmonate (MJ), silver nitrate (SN) and ammonium ceric sulphate (ACS) were also compared for the first time, and the results revealed that expression of TmTXS was all induced by the tested three treatments and the induction effect by MJ was the strongest, implying that TmTXS was high elicitor responsive.